Elliptic flow fluctuations in heavy ion collisions and the perfect fluid hypothesis 
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We analyze the recently measured V2 fluctuation in the context of establishing the degree of fluidity 
of the matter produced in heavy ion collisions. We argue that flow observables within systems with 
a non-negligible mean free path should acquire a "dynamical" fluctuation, due to the random nature 
of each collision between the system's degrees of freedom. Because of this, V2 fluctuations can be 
used to estimate the Knudsen number of the system produced at the relativistic heavy ion collider 
(RHIC). To illustrate this quantitatively, we apply the UrQMD model, with scaled cross sections, 
to show that collisions at RHIC have a Knudsen number at least one order of magnitude below the 
expected value for an interacting hadron gas. Furthermore, we argue that the Knudsen number is 
also bound from below by the V2 fluctuation data, because too small a Knudsen number would break 
the observed scaling of V2 fluctuations due to the onset of turbulent flow. We propose, therefore that 
V2 fluctuation measurements, together with an understanding of the turbulent regime for relativistic 
hydrodynamics, will provide an upper as well as a lower limit for the Knudsen number. We also 
argue that an energy scan of V2 fluctuations could shed light on the onset of the fluid regime. 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq 



One of the most widely cited news (both in the aca- 
demic and popular press) coming out of the heavy ion 
community concerns the discovery of a "perfect fluid" 
in collisions of heavy ions at RHIC 1-5]. The evidence 
for this claim comes from the successful modeling of the 
anisotropic expansion of the matter in the early stage of 
the reaction by means of ideal hydrodynamics [6-8] . This 
argument is compounded by the sensitivity of anisotropic 
expansion to shear viscosity 0, |T(| ■ The presence of a 
non- negligible shear viscosity, therefore, can be detected 
by a careful analysis of anisotropic expansion data. 

However, apart from this evidence for a small viscosity 
many fundamental properties of the fluid are unknown. 
E.g. it is currently discussed whether the observed fluid 
is a strongly interaction Quark-Gluon Plasma (sQGP) 
(Tlj . a bound state Quark-Gluon Plasma (bsQG P) fl3| 
or a (turbulent) Glasma [l3j with instabilities [H], Il5|. 

The connection between theory and experiment rests 
mainly on a single observable for the anisotropic expan- 
sion, namely the elliptic flow coefficient V2 ■ The parame- 
ter V2 is the second Fourier component of the azimuthal 
anisotropy of the particle momenta given by (Tol [T7| 



v 2 = (cos[2(0 - $rp)]> 



(1) 



where <fi denotes the azimuthal angle of one outgoing par- 
ticles and $rp is the azimuthal angle of the reaction 
plane. The angular brackets denote an average over all 
considered particles from all events. 

It should be stressed that (differential) studies of the 
collective flow are among the earliest predicted observ- 
able to probe heated and compressed nuclear matter [l8[ . 
As the transverse flow is intimately connected to the 
pressure gradients in the early stage of the reaction, it 
provides information on the equation of state (EoS) and 
might therefore be used to search for abnormal matter 
states and phase transitions (l9l - [2l| . 



The elliptic flow is of special importance, because it 
is "self-quenching" 0, [l?], H^]: The angular pressure 
gradients creating the anisotropy extinguish themselves 
shortly after the start of the hydrodynamic evolution. 
Thus, the final is insensitive to later stages of the 
evolution, providing a key hole to the hottest, best ther- 
malized, and possibly deconfined phase of the reaction. 

In this paper we use the recently measured event-by- 
event fluctuations of vi [23|, [24| to further investigate the 
properties of the fluid created at RHIC energies. The 
experimental data suggests that the V2 fluctuations follow 
the fluctuations in initial eccentricity e 



UJ,. 



This relation follows in a straight-forward fashion from 
the proportionality between V2 and the eccentricity in- 
ferred from ideal [fUl] and viscous 9| boost-invariant 
hydrodynamics 



v 2 = (3e 



(3) 



where /3 is approximately constant (Indeed, since e is a 
small dimensionless parameter driving anisotropy, this 
relation can be understood simply in terms of Taylor ex- 
pansion). 

As long as we are far away from the turbulent regime 
(on which we comment later), the deterministic nature of 
hydrodynamics, its applicability event by event and Eq. 
|3] constrain [25[ the effect hydrodynamic propagation has 
on initial state fluctuations to the form 



Sv 2 



initial 



(4) 



Explicit calculations have confirmed that this is the case 
for ideal hydrodynamics [26] . A non-zero viscosity should 
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not alter the proportionality, but just lower the value of 

It is, however, surprising that initial conditions be the 
only source of fluctuations. If the system is treated as a 
collection of interacting particles, the random nature of 
each interaction should add a dynamical component to 
the fluctuation of any flow variable, which depends not 
on the initial conditions (with which it is not correlated) 
but on the random nature of each microscopic collision 

((5v 2 f)cf((6ef)+A% n (5) 

One can quantify the deg ree of perfection of the fluid by 
the Knudsen number [25| defined as: 

T AJ 

^ Jv collisions 

with A the mean free path of the particles, L the typical 
length scale of the system N part i C ies the total number 
of particles and N co m S i ons denoting the total number of 
interactions (soft and hard). 

If the Knudsen number is zero, the system becomes 
a perfect fluid. In this case, flow observables are fully 
deterministic. Hence, the probability of a i>2 at a certain 
time assuming a given eccentricity is a 5 — function. 

P{v 2 \e)~8{v 2 -Pe) (7) 

and hence Ad yn is zero. 

In the limit of large numbers of collisions, correlations 
between collisions become weak, so the probability dis- 
tribution in Eq. [7] becomes 

P {v 2 \e. Kn) ~ _L e -(«2-^) 2 /2<x 2 / 8) 

27T<7 

where a(Kn) goes to as the Knudsen number goes to 
zero. Thus, it is sensible to Taylor expand around Kn, 
so 

A dyn ~a^Kn^ + 0{Kn 2 ) (9) 

Note that the only "small parameter" here is the Knud- 
sen number. All dependence on the nature of degrees 
of freedom and their interactions (in particular, whether 
the particles interacting are hadrons or partons, what is 
the equation of state etc.) is encoded within the parame- 
ter a, which by naturalness is of order unity. In the case 
of a vanishing Knudsen number (the ideal hydrodynamic 
limit) Adyn should vanish. 

Such a scaling, apparent in Kinetic theory, can be also 
derived within hydrodynamics [271 ]: Fluctuations in flu- 
ids include a thermal fluctuation term (irrelevant here 
since i>2 is defined in a way that makes it independent 
of random multiplicity fluctuations) and a dynamical au- 
tocorrelation of the energy-momentum tensor. This au- 
tocorrelation scales linearly with the shear and bulk vis- 
cosity, which in turn depend linearly on the mean free 
path [27[ and the inverse of the typical number of col- 
lisions per particle. Subsequent developments [H, [2j| 



have not altered these basic conclusions, which have also 
been compared to Boltzmann equation simulations (30j . 
The latter comparison makes us confident of the "univer- 
sality" of our sca ling , since hydrodynamic fluctuations 
treated in [27], H^, Ho] concern systems where the Boltz- 
mann equation fails (e.g. water). 

This is important, since the definition in Eq. [5] looks 
more natural within a Boltzmann equation formulation, 
which is in turn based on the scattering approximation 
between interactions. It is not clear whether this is a 
good approximation to use within RHIC. While some 
groups have managed to bring models based on these as- 
sumptions in agreement with RHIC data [3lj , the appear- 
ance of fields and off-shell effects in the strongly coupled 
limit is not unreasonable [32] . This matter is complicated 
by the fact that everyday fluids^ (where compressibility is 
typically large and correlated with viscosity due to Pauli 
blocking effects) are fundamentally different from ultra- 
relativistic ones (where the Pauli principle is not thought 
to be relevant, compressibility is bound by causality and 
typically unrelated to viscosity). In the latter, the every- 
day definition of fluid vs gas (based on compressibility) 
becomes inapplicable, and a small viscosity alone is no 
guarantee that the Boltzmann equation approximation is 
not a good one. 

Our definition of Knudsen number, however, is general 
enough to be independent of these considerations: The 
Knudsen number is simply the ratio of a microscopic scale 
(where quantum randomness is important) to the macro- 
scopic scale, in this case the total size of the system. It 
is easy to see that, to leading order, the ratio of the two 
quantities must generally be ~ n/ (sTL) ~ Kn [33]. Such 
a definition of the Knudsen number allows us to recover 
the scaling found in [27i - [29l j (derived with systems where 
the Boltzmann equation is not a good approximation, 
such as water, in mind). 

It is therefore apparent that u> V2 is a test for the hy- 
pothesis that the system at RHIC is a "perfect fluid" ,i.e. 
a locally thermalized system, where "many" particles un- 
dergo "many" collisions over a "small" fraction of the 
system's evolution. Potentially, this test is considerably 
more model-independent than a hydrodynamic analysis 
of (1*2), since Ad yn scales directly with the Knudsen num- 
ber, and all other factors are of order unity. 

Deviations from this limit, including plasma instabil- 
ities, or clustering, should therefore contribute fluctua- 
tions to v-i [34[ that can be probed by comparison to 
the newly available experimental data. Interestingly, 
the "opposite" limit to hydrodynamics, a classical non- 
Abelian field such as the "color glass condensate" (CGC) 
[3^ . is also fully deterministic, and hence would exhibit 
Adyn = 0. Just like ideal hydrodynamics, however, the 



The usual linguistic usage is "liquid" for the definition based on 
compressibility, and "fluid" for the definition based on viscosity. 
These two are often used as synonym in everyday English. 
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absence of dynamical fluctuations in the CGC is an arti- 
fact of it being an effective description with a zero "small 
parameter" . Here, the "small parameter" giving rise to 
fluctuations would be the inverse of the occupancy num- 
ber of each quantum state. As the initial occupancy num- 
ber at RHIC is 3-4 [H, diminishing to < 1 as the CGC 
melts, we expect large dynamical fluctuations also in a 
CGC prethermal stage, through calculating them is best 
left to a future work. 

In this work, we quantitatively assess the sensitivity of 
u> V2 on Kn by a string Aiadron transport approach, the 
UrQMD v2.3 model H3- To explore the different 
regimes, we rescale the total interaction cross sections 
by factors of 1/2, 1 and 3 to vary the strength of the 
interaction. 

Note that we are not using this rescaled UrQMD as a 
realistic model of the system, but rather as a "toy model" 
to study the scaling of i>2 fluctuations of the Knudsen 
number. We believe this is an appropriate approach for 
the present study because the results can, in a straight- 
forward fashion, be converted into an estimate for the 
Knudsen number in heavy ion collision at RHIC. We ex- 
pect that an analysis with partonic degrees of freedom 
will yield the same scaling with the Knudsen number, 
and a quantitative result within the same order of mag- 
nitude, the differences being encoded in the constant a 
of Eq[9] 

As UrQMD is a quantum molecular dynamics simula- 
tion, the Knudsen number can be effectively "measured" 
by keeping track of the collisions during the system's life- 
time, with Kn — {N par ticies /N co ui S i ons ) — 0.6-1.5. This 
analysis parametrically agrees with an estimate following 
(25j . from the ratio of the calculated elliptic flow to the 
hydrodynamic expectation. 

UrQMD also accounts for the expected non-flow ef- 
fects, as well as the fluctuations in the initial condition, 
that also contribute to u V2 [38] . For a general discus- 
sion of the V2 analysis within this approach the reader is 
referred to [H, H3| ■ 

The results of the present calculations are shown in 
Fig. [1] As can be seen, lj V2 and Kn have the expected 
qualitative dependence on the over-all scaling parameter: 
As the factor used to rescale interaction cross sections in- 
creases, (v?) increases [38j and uj V2 decreases. However, 
both (V2) [38[ and w„ 2 are well away from the data-points 
even if the cross section is increased by a factor of three. 
Beyond the given increase, we run into technical diffi- 
culties and grossly over-estimate the total multiplicity of 
the system. Hence, we are not able to explore the scaling 
of the cross-section further than three times the physical 
one within the present approach. 

Fig. [2] shows the scaling of oj V2 w.r.t. the Knudsen 
number. The full line shows a fit of the calculations as- 
suming a the additional fluctuations can be modeled by 
Eqs. [5] and [51 with a and (3 are extracted from the fit. 
While a varies with the fundamental properties of the 
system even in the Poissonian limit, the fitted value is 
sufficient for an order of magnitude estimate of Kn^ 1 . 
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FIG. 1. (color online) UrQMD result for uj v . 2 (panel (a) ) and 
the Knudsen number (panel (b)) as a function of the number 
of participants for Au+Au reaction at \fs NN = 200 GeV. The 
data are taken from [23, l2~ij. 



The scaling in Fig. [2] can be used to extrapolate to 
the inverse Knudsen number needed to describe the el- 
liptic flow fluctuations observed at RHIC. Extrapolating 
the dependence of Eq. [9] to the upper experimental error 
bar for ui V 2 yields a lower bound on the inverse Knudsen 
number on the order of a hundred. This estimate for the 
Knudsen number (or the potential scaling factor for the 
cross section, as the Knudsen number is inversely propor- 
tional to the cross section) qualitatively agrees, with the 
opacity estimate derived from (z^) using pQCD transport 
calculations fill 42| , as well as with the estimate obtained 
through transverse momentum fluctuations [43j . 

With the present calculation, we have established an 
upper limit to the Knudsen number. It is natural to ask 
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what happens to uj V2 as the Knudsen number goes to 
zero, and the system is more closely approximated by 
ideal hydrodynamics. 

The Knudsen number is related to another well known 
number in hydrodynamics, namely Reynold's number. 
Reynold's number is defined as 

tfe=^ = £^M~3# , (10) 
77 77 Kn 

with H being the enthalpy, n being the viscosity, s being 
the entropy density, T denoting the temperature and (v) 
the typical flow velocity. Thus, a small Knudsen number 
goes hand in hand with an increase of Reynold's number. 

However, too high Reynolds numbers inevitably lead 
to instabilities of the hydrodynamic flow (the turbulent 
regime) and will add an additional source of fluctuations 
to uj V2 , due to instabilities in the flow formation. Es- 
timating the Reynold's number for the present trans- 
port simulations leads to Re ~ 1. However using the 
presently advocated ADS/CFT bound r}/s = 1/4tt 0] 
leads to Reynolds numbers well into the ~ 10 2 in the ini- 
tial stages of the hydrodynamic evolution (Fig. for 
T = 200 MeV, L = 10/m and (v) ~ l/%/3 (the speed of 
sound for a relativistic ideal gas), we have Re ~ 100,. 

Following [27j hydrodynamic instabilities will be 
present starting from Re = 10 - 100. If Re > 100 - 1000 
the flow will generally become turbulent [27j . although 
the onset of turbulence will also depend on the bound- 
ary conditions: the larger the bluffness of a layer of fluid 

(defined by ( u l x dA ) , where u l is the flow vector and 



10' 



dA the layer surface element) , the less Reynolds number 
is required for the onset of turbulence. For a compress- 
ible fluid expanding from an "almond-like" shock, the 
last condition is likely to be satisfied close to the "edges" 
of the almond, provided that compressibility does not 
quench the onset of turbulence (the last question is not 
conclusively settled, through recent evidence [45l.l46| sug- 
gests that adding compressibility does not significantly 
change the critical Reynolds number for the onset of tur- 
bulence) . 

Thus, we are led to conclude that, below a certain 
critical Knudsen number, the scaling in Eq. [5] should 
break down and w„ 2 should increase significantly above 
the "ideal" Eq. Rvalue: While Ad yn would continue to 
decrease with decreasing Knudsen number, uj V2 would not 
anymore scale with ((5e) 2 ) / (e) but the initial fluctuation 
would be amplified by the turbulent evolution. If t is the 
lifetime of the system and tq the timescale of the evolu- 



tion of turbulence, 



in a turbulent fluid would scale 



as an exponent of an approximately power-law function 
of the initial volume (i.e. the number of participants) 



OJ,, 



<(M 2 ) 



"2 



,t/t 



<(^) 2 > 



exp [N } 



part] 



(11) 



From causality (the time it takes for rarefaction waves to 

travel across the system, ~ size/c s ~ Np^ t /c s , where c s 
is the speed of sound), it can be deduced that n ~ 1/3. 
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FIG. 2. (color online) Relationship between uj V2 and the 
Knudsen number, plotted together with the Poissonian ex- 
pectation. The parameter a was fitted from the data. See 
Fig. [I] for the legend 



It is apparent from Fig. Q]that such a scaling is not ob- 
served in the experimental data [23|, [24| , so the viscosity 
of the system created at RHIC is high enough to place it 
out of the turbulent regime. This sustains the argument 
that the mere observation of a well-defined (1)2) places a 
lower constraint on viscosity because it signals that the 
system is not in a sufficiently turbulent regime. 

The implications of this statement on the closeness of 
the fluid created at RHIC to the ADS/CFT viscosity 
bound are still not clear. It is difficult to make a more 
precise estimate since the turbulence in the system pro- 
duced in heavy ion collisions has not as yet been studied 
(for first attempts with QCD transport approaches, the 
reader is referred to [H, [ll, S3, El). 

In one dimension, the stability of boost-invariant dy- 
namics (the boundary condition used, either exactly or in 
approximate form for simulations at RHIC energy) has 
been thoroughly studied. Boost-invariant evolution was 
found to be generally stable at the early stages, where vi 
forms 0, [Hc| (thr oug h instabilities could play a big role 
during freeze-out [5l[ |52]|). This leads us to think that 
if the system does have an early turbulent stage, tq re- 
mains long compared to its duration. This conclusion is 
however bound to change within full 3D hydrodynamics, 
especially if the system is not to a good approximation 
boost invariant, as recent initial state calculations sug- 
gest [nm. 

Thus, before a quantitative answer to these questions 
can be given, a transport or hydrodynamic model capable 
of modeling turbulence at the scale of heavy ion collisions, 
and hence of inferring a quantitative value of tq in Eq. 
111! is necessary. Up to now, the only known calculation 
of the Reynolds number and the onset of turbulent flow 
in heavy ion collision has been done in Ref. [oH ]. 

The onset of turbulence could be signaled experimen- 
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FIG. 3. (color online) The Reynolds number of a system the 
size of a collision between nuclei A, with the viscosity given 
by the conjectured ADS/CFT universal bound. 



tally by a widening of ui V 2 and a change of its dependence 
on Np ar t from constant to exponential scaling as per Eq. 
[TT1 Thus, combined with the data in the non-turbulent 
regime, analyzed using the ansatz of Eq. O the experi- 
mental measurement of lj V2 in a wide range of energies 
and system sizes can yield a lower as well as an upper 
limit of K n. 

Moreover, the energy dependence of uj V2 could acquire 
a crucial phenomenological role in the light of the uni- 
versal scaling seen in v 2 /e [54l - [5ti| (panel (a) in Fig. @|. 
The scaling variable is the multiplicity rapidity density 
normalized by the initial overlap surface dN / [Sdy), cho- 
sen because it corresponds, in the boost-invariant picture 
(HI, to the entropy density divided by thermalization 
time. One could interpret this scaling as the approach 
to the ideal hydrodynamics limit as the initial density 
become large. If this interpretation is correct, the Knud- 
sen number smoothly decreases as an inverse power of 
dN/(Sdy), but has little sensitivity to the change in de- 
grees of freedom at the phase transition [57) . Alterna- 
tively, it could be that the large flow observed at RHIC is 
indicative of a downward "jump" in the Knudsen number 
when the critical initial density needed to free partonic 
degrees of freedom is achieved. 

The observation of lu V2 , and the excitation function of 
uj V2 /lu <l could differentiate between these scenarios (Fig. 
0] panel (b)). If the system smoothly becomes more fluid 
at greater density, u> V2 /uj e can be expected to decrease in- 
versely with dN / (Sdy) (smoothly if this is a continuous 
approach to hydrodynamics or abruptly if the a transi- 
tion to a more fluid regime is linked to a phase transi- 
tion). If fluidity is present in all systems to the same 
amount, uj V2 /uj e will stay constant across energies and 
system sizes. 

In conclusion, we have argued that the experimental 




Transition: fluid phase 
Hydrodynamic regime always there 
dN/(Sdy) 

FIG. 4. (color online) A qualitative plot (panel (b)) show- 
ing the scaling of the V2 fluctuations in the three scenarios 
suggested by the scaling of Vi across energy and system size 
[5j-[53], shown in panel (a). The x-axis, the rapidity density 
normalized by the overlap area, corresponds to the entropy 
density in the Bjorken hydrodynamic scenario. If the hydro- 
dynamic limit is smoothly approached with increasing system 
volume/lifetime the difference between the observed t>2 fluctu- 
ation and the initial estimate should go as the red dashed line. 
If the hydrodynamic limit indicates a transition between a vis- 
cous hadronic gas and the sQGP, the scaling with 1/ SdN/dy 
should be broken, with higher energy ("sQGP regime") lower 
centrality events having a lower V2 fluctuation than equivalent 
more central ("hadronic regime") events. This is indicated in 
the plot by the blue dot-dashed line. Finally, a constancy of 
uj V 2 /us t might indicate that the "hydro regime" was actually 
with us all along, and only the initial conditions are respon- 
sible for the apparent rise in V2- 



observation of u> V2 can provide unique information to esti- 
mate the Knudsen number Kn, and hence to to pin down 
the perfection of the fluid created in heavy ion collisions 
quantitatively. We have used a transport model to esti- 
mate a lower limit of Kn~ l , and found that it is nearly 
two orders of magnitude below the value needed to de- 
scribe the v% fluctuations at RHIC. We have also argued 
that the currently observed scaling of lu V2 should break 
in the turbulent regime, and hence the measurement of 
w„ 2 potentially places an upper as well as a lower limit 
on Kn~ x . We have furthermore suggested that, in the 
light of these considerations, an energy and system size 
scan of f 2 fluctuations can shed light on the approach to 
the hydrodynamic regime. However, before these limits 
can be quantitatively ascertained, much more theoretical 
modeling and experimental investigation is required. 
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